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Abstract 

We report an improved measurement of direct CP violation in the decay K~^Tr~, and 

a search for CP violation in the decays — > ii'+Tr^ and -6+ — -k'^tt^. The measured CP 
violating asymmetries are AcpiK'^'n-~) = -0.113 ± 0.022(stat.) ± 0.008(syst.), .4cp(^+7r°) = 
0.04 ± 0.04(stat.) ± 0.02(syst.) and ^cp(7r+7r°) = 0.02 ± 0.08(stat.) ± 0.01(syst.), where the latter 
correspond to the intervals -0.03 < Acp{K^Tr°) < 0.11 and -0.12 < ^C7p(vr+vr°) < 0.15 at 90% 
confidence level. These results arc obtained from a data sample that contains 386 million BB pairs 
that was collected near the T(45') resonance, with the Belle detector at the KEKB asymmetric 
energy e'^e~ collider. All of the results are preliminary. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh, 14.40.Nd 
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In the Standard Model (SM) CP violation arises via the interference of at least two 
diagrams with comparable amplitudes but different CP conserving and violating phases. 
Mixing induced CP violation in the B sector has been established in 6 — ccs transitions 
P, 01 . Direct CP violation is expected to be sizeable in the B meson system j^. The 
first experimental evidence for direct CP violation was shown by Belle for the decay mode 
B^ — > vr+vr" 01 • This result suggests large interference between tree and penguin diagrams 
and the existence of final state interactions In addition, both Belle Q and BABAR 0] 
have recently reported evidence for direct CP violation in the decay B^ — K^t:^ . 

The partial rate CP violating asymmetry is defined as: 

^ iV(g ^ 7) - N{B ^ /) 
'''' N{B^f) + N{B^fy 

where N{B —>■ f) is the yield for the B Kti/titi decay and N{B /) denotes that of 
the charge-conjugate mode. Theoretical predictions with different approaches suggest that 
Acp{K~^'7i~) could be either positive or negative jSj]. Although there are large uncertain- 
ties related to hadronic effects in the theoretical predictions, results for AcpU^^'n'^) and 
AcpiK~^n^) are expected to have the same sign and comparable magnitudes j8|. However, 
our previous measurements show that Acp{K~^7i~) and Acp{K~^7i^) are opposite in sign 
(although AcpiK+w'^) is consistent with no asymmetry), and their central values are found 
to deviate from each other by 2.4cr. These findings are consistent with those reported by 
BABAR 0, 1^ • It is suggested that the disagreement may be due to the contribution of the 
electroweak penguin process or other mechanisms Therefore, it is important to verify 
whether the discrepancy persists with improved precision. In this Letter, we report Acp 
measurements using 386 million BB pairs collected with the Belle detector at the KEKB 
e~^e~ asymmetric-energy (3.5 on 8 GeV) collider jll| operating at the T(4S') resonance. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electro- magnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere fl5]. Two inner detector configurations 
were used. For the first sample of 152 million BB pairs (Set I), a 2.0 cm radius beampipe 
and a 3-layer silicon vertex detector were used; for the latter 234 million BB pairs (Set II), a 
1.5 cm radius beampipe, a 4- layer silicon detector and a small-cell inner drift chamber were 
usedlli 



The B candidate selection is the same as that described in Ref. [1^. Charged tracks are 
required to originate from the interaction point (IP). Charged kaons and pions are identified 
using dE/dx information and Cherenkov light yields in the ACC. The dE/dx and ACC 
information are combined to form a K-tt likelihood ratio, TZ{Ktt) = Ck/{Ck + ^tt), where 
Ck/tt is the likelihood of kaon/pion. Charged tracks with TZ^Kn) > 0.6 are regarded as 
kaons and TZ{Ktt) < 0.4 as pions. Furthermore, charged tracks that are positively identified 
as electrons are rejected. The K/n identification efficiencies and misidentification rates are 
determined from a sample of kinematically identified D*^ — * D^tt^, — > K~tt'^ decays. 
Table m shows the results. It is clear that the detection efficiency for K~ /n^ is greater than 
for K~^/'ir~] this small efficiency bias will be corrected in the Acp measurement. 

Candidate tt" mesons are reconstructed by combining two photons with invariant mass 
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TABLE I: Performance oi K — it identification measured using D*^ D^ir^, K vr^ decays. 



Set I Set II 

Eff. (%) Fake rate (%) Eff. (%) Fake rate 

if+ 83.76 ±0.18 5.10 ±0.12 81.92 ±0.15 6.29 ±0.10 

ET" 84.76 ±0.18 5.69 ±0.12 82.79 ±0.14 6.71 ±0.10 

7r+ 91.24 ±0.15 10.72 ±0.15 89.88 ±0.12 12.28 ±0.12 

TT" 90.53 ±0.15 10.08 ±0.15 89.08 ±0.13 11.83 ±0.12 



between 115 MeV/c^ and 152 MeV/c^, which corresponds to ±2.5 standard deviations. Each 
photon is required to have a minimum energy of 50 MeV in the barrel region (32° < 9^ < 
129°) or 100 MeV in the end-cap region (17° < 0^ < 32° or 129° < 0^ < 150°), where 9^ 
denotes the polar angle of the photon with respect to the beam line. To further reduce 
the combinatorial background, 7r° candidates with small decay angles (cos^* > 0.95) are 
rejected, where 9* is the angle between a vr'' boost direction from the laboratory frame and 
its 7 daughters in the 7r° rest frame. 

Two variables are used to identify B candidates: the beam-constrained mass, Mbc = 

\/-^beam ~ 'P*B i ^^^^ energy difference, /S.E = — -Ebeam; where -Ebeam the beam 
energy and E'^ and p*^ are the reconstructed energy and momentum of the B candidates in 
the center-of-mass (CM) frame. Events with Mbc > 5.20 GeV/c^ and -0.3 GeV < < 
0.5 GeV are selected for the final analysis. 

The dominant background is from e~^e~ qq {q = u, d, s, c) continuum events. To 
distinguish the signal from the jet-like continuum background, event topology variables and 
B flavor tagg ing information are employed. We combine a set of modified Fox- Wolfram 
moments [l5| into a Fisher discriminant. The probability density functions (PDF) for this 
discriminant, and the cosine of the angle between the B flight direction and the z axis, are 
obtained using signal and continuum Monte Carlo (MC) events. These two variables are 
then combined to form a likelihood ratio 71 = Cs/{Cs + jCgg), where Cs(qq) is the product 
of signal {qq) probability densities. Additional background discrimination is provided by 
B flavor tagging. The standard Belle flavor tagging algorithm gives two outputs: a 
discrete variable indicating the flavor of the tagging B and the MC-determined dilution 
factor r, which ranges from zero for no flavor information to unity for unambiguous flavor 
assignment. An event that contains a lepton (r close to unity) is more likely to be a BB 
event so a looser TZ requirement can be applied. We divide the data into r > 0.5 and r < 0.5 
regions. The continuum background is reduced by applying a selection requirement on TZ 
for events in each r region of Set I and Set II according to the figure-of-merit defined as 
N^^P/ ^ Ns^^ + Ng^^, where A^^^^ denotes the expected signal yields based on our previous 
branching fraction measurements and A^|^ denotes the expected qq yields from sideband 
data (Mbc < 5.26 GeV/c^). A typical requirement suppresses 92-99% of the continuum 
background while retaining 48-67% of the signal. 

Backgrounds from T{AS) —>■ BB events are investigated using a large MC sample. Af- 
ter the TZ requirement, we find a small charmless three-body background at low AE, and 
reflections from other B^ vr+vr" decays due to K-n misidentification. 

The signal yields are extracted by applying unbinned two dimensional maximum likeli- 
hood (ML) fits to the (Mbc, AE) distributions of the B and B samples. The likelihood for 
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each mode is defined as 



L = exp (- Y.^s,Ki) X I[iT. Ns,k,J'Ps,k,J,^) (2) 

Vs,k,j,i = - (li'AcPj\Ps,k,o{Mhci. AEi), (3) 

where s indicates Set I or Set II, k distinguishes events in the r < 0.5 or r > 0.5 regions, i 
is the identifier of the i-th event, P(Mbc, Ai?) is the two-dimensional PDF, q indicates the 
reconstructed B meson fiavor, B{q = +1) or B{q = — 1), Nj is the number of events for the 
category j, which, in turn, corresponds to either signal, qq continuum, a refiection due to 
K-TT misidentification, or background from other charmless three-body B decays. 

The yields and asymmetries for the signal and backgrounds are allowed to fioat in all 
modes. Since the i^+vr" and tt+tt" refiections are difficult to distinguish with AE and Mbc 
we fit these two modes simultaneously with a fixed refiection-to-signal ratio based on the 
measured K-n identification efficiencies and fake rates. All the signal PDFs (P(Mbc, A£')) 
are obtained using MC simulations based on the Set I and Set II detector configurations. 
The same signal PDFs are used for events in two different r regions. No strong correlations 
between Mbc and AE are found for the B K^n^ signal. Therefore, its PDF is modeled 
by a product of a single Gaussian for M^c and a double Gaussians for AE. For the modes 
with a 7r° meson in the final state, there are correlations between M^c and AE in the tails 
of the signals; hence, their PDFs are described by smoothed two-dimensional histograms. 
Discrepancies between the signal peak positions in data and MC are calibrated using B~^ —>■ 
£)°7r+ decays, where the D° — > K'^n'n^ sub-decay is used for the modes with a 7r° meson 
while D° — * K~^TT~ is used for the K~^7r~ mode. The MC-predicted AE resolutions are 
verified using the invariant mass distributions of high momentum D mesons. The decay 
mode — > K~^tt~ is used for B^ — > K~^7t~, and — > K+7r~7r° for the modes with a 7r° in 
the final state. The parameters that describe the shape of the signal PDFs are fixed in all 
of the fits. 

The continuum background in AE is described by a first or second order polynomial while 
the Mbc distribution is parameterized by an Argus function f{x) = x\/l — x'^ exp [—^(1 — 
x^)], where x is Mbc divided by half of the total center of mass energy. The continuum 
PDF is the product of an Argus function and a polynomial, where parameters ^ and the 
coefficients of the polynomial are free parameters. These free parameters are r-dependent. 
A large MC sample is used to investigate the background from charmless B decays and 
a smoothed two-dimensional histogram is taken as the PDF. The functional forms of the 
PDFs are the same for the B and B samples. 

The efficiency of particle identification is slightly different for positively and negatively 
charged particles; consequently, the parameter Acp in Eq.3 becomes an effective partial rate 
asymmetry. For the K'^tt^ and tt+tt^ modes, this raw asymmetry can be expressed as: 

-^^^ - 1 + A.Acp' ^ ' 

where Acp is the true partial rate asymmetry and the efficiency asymmetry Ae is the 
efficiency difference between K~{'it~^) and K~^{7r~) divided by the sum of their efficiency. 
The situation is more complicated for the K~^7r~ mode because, in addition to the bias due 
to the efficiency difference between K~n~^ and K~^7i~ , a K~7t~^ signal event can be doubly 
misidentified as a K'^7r~ candidate and dilute Acp- The efficiency asymmetry results in a 



6 




AE (GeV) AE (GeV) 



FIG. 1: Mbc (top) and AS (bottom) distributions for B K'tt+ (left) and K+tt~ (right) 

candidates. The histograms represent the data, while the curves represent the various components 
from the fit: signal (dashed), continuum (dotted), three-body B decays (dash-dotted), background 
from misidentification (hatched), and sum of all components (solid). 

Acp bias of +0.01, while the small dilution factor due to double misidentification reduces 
the Acp by a factor of 0.98. 

Table ITTl shows the signal yields and Acp values for each mode. The asymmetries for the 
background components are consistent with zero within errors. Projections of the fits are 
shown in Figs. 1-3. The systematic errors from fitting are estimated from the deviations in 
Acp after varying each parameter of the signal PDFs by 1 standard deviation. The uncer- 
tainty in modeling the three-body background is studied by excluding the low AE region 
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FIG. 2: Mbc (top) and t^E (bottom) distributions for B' R-tt^ (left) and 5^ 
candidates. The curves are described in the caption of Fig^ 



K+TT^ (right) 



(< —0.12 GeV) and repeating the fit. Systematic uncertainty due to particle identification 
is estimated by repeating the fit after varying the K/ tt efficiencies and fake rates by 1 stan- 
dard deviation. At each step, the deviation in Acp is added in quadrature to provide the 
systematic errors, which are less than 0.01 for all modes. A possible bias from the fitting 
procedure is checked in MC and a bias due to the TZ requirement is investigated using the 
— > D^TT^ samples. No significant bias is observed. The systematic uncertainties due 
to the detector bias are tested using the fit results for the continuum background listed in 
Table |n] We find a small background asymmetry dependence on the TZ requirement for the 
K^T[~ mode, and assign the uncertainty from the fit result of the D^7c~^{D^ —>■ K^ti~) 
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TABLE II: Fitted signal yields, Acp results, and background asymmetries for individual modes. 



Mode 


Signal Yield 


Acp 


Bkg Acp 




3026 ± 63 


-0.113 ±0.022 ±0.008 


-0.001 ± 0.004 




1084 ± 45 


±0.04 ±0.04 ±0.02 


-0.02 ±0.01 




454 ± 36 


±0.02 ±0.08 ±0.01 


-0.01 ±0.01 



sample (±0.007) as the systematic uncertainty due to detector bias. The final systematic 
errors are then obtained by quadratically summing the errors due to the detector bias and 
the fitting systematics. 

The partial rate asymmetry Acp{K^tt^) is found to be —0.113 ± 0.022 ± 0.008. The 
significance including the effect of systematic uncertainty is A.97a. This result supersedes 
our previous measurement 0] and remains consistent with the value reported by BABAR, 
AcpiK+n-) = -0.133 ± 0.030 ± 0.009 0. The observed Acp{K+n'^) value is consistent 
with zero at the current level of statistical precision. The difference between the results for 
Acp{K^n^) and AcpiK^n'^) persists; their central values differ by S.lcr. This suggests a 
possible contribution from the electroweak penguin process or other mechanisms 3|. No 



evidence of direct CP violation is observed in the decay vr+vr^. We set 90% CL. 

intervals: -0.03 < Acp{K+n'^) < 0.11 and -0.12 < Acpin+n'^) < 0.15. All of the above 
results are preliminary. 
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